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ABSTRACT:  We  analyzed  the  thermodynamic  and  structural 
determinants  of  indolicidin  interactions  with  eukaryotic  and 
prokaryotic  cell  membranes  using  a  series  of  atomistically 
detailed  molecular  dynamics  simulations.  We  used  quartz- 
supported  bilayers  with  two  different  compositions  of 
zwitterionic  and  anionic  phospholipids  as  model  eukaryotic 
and  prokaryotic  cell  membranes.  Indolicidin  was  preferentially 
attracted  to  the  model  prokaryotic  cell  membrane  in  contrast 
to  the  weak  adsorption  on  the  eukaryotic  membrane.  The 
nature  of  the  indolicidin  surface  adsorption  depended  on  an 
electrostatic  guiding  component,  an  attractive  enthalpic 
component  derived  from  van  der  Waals  interactions,  and  a 
balance  between  entropic  factors  related  to  peptide  confine¬ 
ment  at  the  interface  and  counterion  release  from  the  bilayer 
surface.  Thus,  whereas  we  attributed  the  specificity  of  the 
indolicidin/membrane  interaction  to  electrostatics,  these 
interactions  were  not  the  sole  contributors  to  the  free  energy  of  adsorption.  Instead,  a  balance  between  an  attractive  van  der 
Waals  enthalpic  component  and  a  repulsive  entropic  component  determined  the  overall  strength  of  indolicidin  adsorption. 


■  INTRODUCTION 

Antimicrobial  peptides  are  small  cationic  peptides  that  play 
important  roles  in  the  host-defense  mechanism  of  many  living 
organisms  against  bacterial  infections.1-3  Antimicrobial  pep¬ 
tides  act  rapidly  against  a  broad  spectrum  of  microbes  and  have 
attracted  attention  as  potential  novel  antibiotics  to  combat 
multidrug-resistant  bacterial  infections,  either  as  standalone 
agents  or  in  combination  therapies  with  other  conventional 
antibiotics.4-6  Indolicidin  is  an  antimicrobial  peptide  isolated 
from  bovine  neutrophils.7  It  consists  of  13  amino  acid  residues 
(NH3+-ILPWKWPWWPWRR-NH2)  with  an  uncapped  N- 
terminus  and  an  amidated  C -terminus  resulting  in  a  net 
positive  charge  of  4  e.  The  structure  of  indolicidin  is  known  to 
be  disordered  in  aqueous  solution  but  has  been  shown  to  form 
an  extended  conformation  when  bound  to  detergent  micelles.8 
The  mechanism  of  indolicidin  action  against  both  bacterial  and 
human  cells  has  been  investigated  both  experimentally9-17  and 
computationally.18-20  Even  though  several  different  mecha¬ 
nisms  of  action  by  indolicidin  have  been  proposed,  the  exact 
mechanism  remains  unclear.14,20  The  delineation  of  the 
mechanism  is  not  straightforward  because  the  biological  action 
that  ultimately  results  in  lysis  of  bacteria  is  multifaceted. 
Initially,  indolicidin  must  recognize  and  attach  to  bacterial  cell 
membranes.  After  this  attachment,  several  lytic  mechanisms  are 
possible,  for  example,  membrane  thinning  resulting  in  a 
structural  weakening  of  the  membrane  and  loss  of  cell 
homeostasis,  interactions  with  proteins  in  the  membrane,  or 
translocation  to  the  cytosol  where  the  peptide  can  further 
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interact  with  diverse  cellular  targets,  such  as  RNA  and  DNA,  or 
other  highly  charged  species.  Existing  results  from  both 
experimental8,14,16  and  computer  simulation18-20  work  confirm 
that  the  indolicidhfs  mechanism  of  action  includes  an 
important  localization  component  associated  with  peptide 
adsorption  to  the  solvent  side  of  membrane  bilayer  interface. 
In  this  work,  we  quantify  and  address  the  detailed 
thermodynamic  and  structural  elements  of  this  initial  attraction 
of  indolicidin  to  model  prokaryotic  and  eukaryotic  bilayers 
using  computational  approaches. 

Solid-supported  lipid  bilayers  are  widely  used  as  experimental 
model  systems  to  investigate  detailed  interactions  of  biological 
molecules  with  membranes.  Such  model  systems  have  been 
used  to  investigate  the  interactions  of  antimicrobial  peptides 
with  lipid  membranes  using  quartz  crystal  microbalance 
(QCM)  measurements,12,21  sum  frequency  generation  (SFG) 
vibrational  spectroscopy,22,23  and  in  situ  atomic  force 
microscopy  (AFM). 14-16  Similarly,  one  can  use  computational 
models  to  create  realistic  solid-supported  lipid  bilayers  at  the 
atomic  level24-27  to  map  out  the  interactions  of  antimicrobial 
peptides  with  lipid  membranes  in  great  detail.  Recently,  we 
successfully  modeled  solid-supported  bilayers  with  a  hydrated 
dimyristoylphosphatidylcholine  (DMPC)  lipid  bilayer  depos¬ 
ited  on  a  quartz  crystal  surface.28  Importantly,  we  showed  a 
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critical  need  to  implement  the  appropriate  electrostatic 
boundary  conditions  for  simulations  and  analyses  to  capture 
the  structural  and  electrostatic  properties  of  DMPC  lipid 
bilayers  that  correspond  to  the  physical  systems  being  modeled, 
for  example,  solvated  bilayers  representing  lamellar  systems, 
nonlamellar  bilayers  mimicking  membranes,  and  quartz- 
supported  bilayers.  These  methodological  considerations  are 
critical  because  the  electrostatic  interaction  between  cationic 
antimicrobial  peptides  and  the  anionic  lipids  in  the  bacterial  cell 
membrane  is  a  key  feature  that  governs  the  selective  attraction 
of  antimicrobial  peptides  to  the  bacterial  cell  membrane.3'29  A 
detailed  understanding  of  this  electrostatic  interaction  needs  to 
take  into  consideration  all  the  complex  interactions  among 
explicit  solvent  molecules,  counterions,  and  salt  ions  present 
under  physiological  conditions.  In  particular,  counterions  and 
salts  strongly  diminish  the  electrostatic  interactions  between 
antimicrobial  peptides  and  anionic  lipids,  leading  to  a  more 
nuanced  interpretation  of  the  electrostatic  effect  associated  with 
cationic  antimicrobial  peptides. 

In  this  work,  we  investigated  the  differences  in  the 
interaction  of  the  antimicrobial  peptide  indolicidin  with 
model  bacterial  and  eukaryotic  cell  membranes  by  performing 
atomistically  detailed  molecular  dynamics  (MD)  simulations 
that  included  solvent  molecules,  salt  ions,  and  lipids.  We 
constructed  model  prokaryotic  and  eukaryotic  cell  membranes 
using  mixed  (zwitterionic  and  anionic)  and  neutral  (zwitter- 
ionic)  lipid  bilayers  respectively  deposited  on  a  quartz  crystal 
surface.  To  characterize  the  interaction  of  indolicidin  with  the 
lipid  membranes,  we  calculated  the  potential  of  mean  force 
(PMF)30  or  free  energy  profile  with  respect  to  the  distance  of 
indolicidin  from  the  lipid  membrane  using  a  series  of  MD 
simulations.  We  confirmed  the  existence  of  the  selective 
attractive  interaction  of  indolicidin  with  the  model  bacterial 
cell  membrane,  which  results  in  indolicidin  being  preferentially 
located  at  the  membrane  interface.  We  further  characterized  the 
nature  of  the  underlying  interaction  in  terms  of  structural  and 
thermodynamic  properties.  In  particular,  we  found  that 
placement  of  indolicidin  at  the  membrane  surface  was 
associated  with  a  minimum  PMF  of  —3.5  ±  0.2  kcal/mol  at 
the  prokaryotic  membrane,  compared  to  —0.6  ±  0.6  kcal/mol 
at  the  eukaryotic  membrane.  We  further  used  the  distance- 
dependent  PMF  to  estimate  the  free  energy  of  indolicidin 
adsorption  (A G°)  to  be  —1.94  ±  0.56  kcal/mol  for  the  model 
prokaryotic  cell  membrane  and  —0.05  ±  0.46  kcal/mol  for  the 
eukaryotic  membrane.  The  detailed  thermodynamic  and 
structural  analysis  of  these  systems  provided  support  for  the 
idea  that  the  specificity  of  the  indolicidin/membrane 
interaction  can  be  attributed  to  electrostatics,  but  these 
interactions  were  not  the  sole  contributors  to  the  free  energy 
of  adsorption.  Instead,  a  balance  between  an  attractive  van  der 
Waals  enthalpic  component  and  a  repulsive  entropic 
component  determined  the  overall  strength  of  indolicidin 
adsorption. 

■  METHODS 

Simulations.  MD  simulations  were  performed  with  the 
NAMD  MD  simulation  program31  using  the  CHARMM  27  all 
atom  force  field.32  We  used  DMPC  and  dimyristoyl 
phosphatidylglycerol  (DMPG)  as  model  zwitterionic  and 
anionic  phospholipids,  respectively.  The  TIP3P  model33  was 
used  to  describe  water  molecules.  The  force  field  for  the  quartz 
(Oil)  crystal  was  taken  from  one  developed  by  Lopes  et  al.34  A 
primitive  unit  cell  of  quartz  (Oil)  of  ~15  A  thickness  was 


replicated  in  two  dimensions  to  construct  a  quartz  (Oil)  crystal 
with  a  surface  area  of  36.5  X  37.0  A2  in  the  unit  cell,  as  in  the 
previous  study  by  Lopes  et  al.34  and  in  our  recent  study.28  One 
side  of  the  constructed  crystal  was  covered  with  hydrophilic 
silanols  (Si— OH),  and  silicon  atoms  on  the  other  side  were 
saturated  with  hydrogens  (Si— H).  Crystal  atoms,  except  for  the 
O— H  groups  of  silanols,  were  held  fixed  to  maintain  the  (Oil) 
crystal  geometry.  Short-range  interactions  outside  a  10  A  cutoff 
were  truncated.  Long-range  electrostatic  interactions  were 
calculated  with  the  particle  mesh  Ewald  method35  with  a 
correction  term  for  the  planar  vacuum  boundary  condition,28'36 
referred  to  as  EW3DC,  which  we  implemented  in  NAMD.  The 
EW3DC  correction  term,  combined  with  a  sufficiently  large  box 
length  ( Lz )  in  the  z  direction  normal  to  the  interface,  effectively 
implements  a  2D  periodic  boundary  condition,  which  is 
appropriate  for  solid-supported  lipid  bilayers.  Bonds  involving 
hydrogens  were  constrained  with  the  SHAKE  algorithm.37  We 
used  a  time  step  of  2  fs  for  the  time  integration,  and  the 
temperature  was  maintained  at  310  K  using  Langevin  dynamics 
with  a  damping  coefficient  of  10  ps-1.38 

Preparation  of  Quartz-Supported  Lipid  Bilayers. 
Quartz-supported  lipid  bilayers  were  prepared  as  illustrated  in 
Figure  1.  Hydrated  lipid  bilayers  are  in  contact  with  the 


Figure  1.  Snapshot  of  the  hydrated  DMPC/G  lipid  bilayer  supported 
on  a  quartz  crystal  illustrates  the  placement  of  the  unit  cell  in  the 
supported  bilayer  system.  The  blue  box  represents  the  unit  simulation 
cell  with  an  Lz  of  250  A.  The  simulation  cell  is  repeated  periodically  in 
all  three  directions.  Gray  and  red  lines  represent  DMPC  and  DMPG 
lipids,  respectively.  Gray  and  red  balls  near  the  membrane  interfaces 
denote  phosphorus  atoms  of  headgroups  in  DMPC  and  DMPG  lipids, 
respectively.  Smaller  pink  spheres  represent  oxygen  atoms  of  water 
molecules.  Blue  and  green  spheres  represent  sodium  and  chloride  ions, 
respectively. 


hydrophilic  side  of  the  quartz  crystal.  Both  a  pure  DMPC 
bilayer  and  a  mixed  DMPC/DMPG  (DMPC/G)  bilayer  were 
prepared  as  model  eukaryotic  and  bacterial  cell  membranes, 
respectively.  For  a  starting  point  of  a  pure  DMPC  bilayer,  we 
used  the  last  configuration  from  the  200  ns  simulation  of  the 
DMPC  bilayer  composed  of  44  DMPC  molecules  (22  on  each 
leaflet)  supported  on  the  quartz  surface  described  in  our 
previous  work.28  Similarly,  a  mixed  lipid  bilayer  with  34  DMPC 
and  10  DMPG  molecules  (17  DMPC  and  5  DMPG  on  each 
leaflet)  supported  on  the  quartz  surface  was  prepared  after  a  45 
ns  equilibration  in  aqueous  solution  and  a  50  ns  equilibration 
on  the  quartz  surface.  An  initial  configuration  of  the  mixed 
DMPC/G  lipid  bilayer  was  prepared  with  the  CHARMM- 
GUI.39  In  these  configurations,  the  upper  water  layer  between 
the  lipid  bilayer  and  the  water/vacuum  interface  contained 
about  twice  as  many  water  molecules  as  the  lower  water  layer 
between  the  bilayer  and  the  quartz.  To  increase  the  bulklike 
region  of  water  to  accommodate  indolicidin,  we  doubled  the 
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number  of  water  molecules  in  the  upper  water  layer.  As  a  result, 
the  upper  water  layer  contains  about  four  times  as  many  water 
molecules  as  the  lower  water  layer.  To  observe  the  effects  of 
salt  ions  on  the  properties  of  lipid  bilayers,  we  performed 
simulations  with  and  without  added  NaCl  salt  for  a  pure 
DMPC  bilayer.  For  a  simulation  with  salt,  we  added  2  and  8 
pairs  of  NaCl  to  the  lower  and  upper  water  layers  respectively 
resulting  in  a  salt  concentration  of  about  0.16  M, 
commensurate  with  physiological  salt  concentrations.  For  a 
simulation  with  the  mixed  DMPC/G  bilayer,  10  extra  sodium 
counterions  (5  to  each  water  layer)  were  added  to  maintain  the 
net  charge  neutrality.  The  area  per  lipid  headgroup  is  about  60 
A2,  which  is  close  to  the  experimental  value  for  DMPC.40  The 
box  length  (Lz)  in  the  z  direction  normal  to  the  interface  was 
fixed  at  250  A,  which  was  large  enough  to  create  a  water/ 
vacuum  interface  as  shown  in  Figure  1.  This  allowed  the  system 
to  adjust  to  an  optimum  density  without  the  need  of  constant- 
pressure  simulations.  Simulation  for  each  prepared  quart- 
supported  lipid  bilayer  system  lasted  100  ns.  The  initial  20  ns 
were  discarded  as  equilibration. 

Preparation  of  Indolicidin.  Five  different  initial  config¬ 
urations  of  indolicidin  were  prepared.  The  first  configuration 
was  taken  from  the  structure  of  indolicidin  bound  to  a 
do decylphospho choline  (DPC)  micelle  as  determined  by  NMR 
(PDB  id:  lg89).8  After  a  brief  energy  minimization,  we 
performed  an  MD  simulation  of  indolicidin  in  aqueous  solution 
for  100  ns.  The  other  four  configurations  of  indolicidin  were 
selected  from  conformations  generated  by  the  ECEPPAK 
program.41'42  We  selected  two  conformations  where  NOE 
distances  deviated  less  than  1  A  from  experimental  values  and 
another  two  conformations  where  the  Ca  rmsd  of  the  core 
residues  3—11  was  less  than  1  A  away  from  the  NMR  model 
structure.  We  performed  an  MD  simulation  lasting  30  ns  in 
aqueous  solution  for  each  of  the  four  selected  conformations. 
We  kept  the  last  conformation  of  indolicidin  from  each  solution 
simulation  and  placed  it  in  the  middle  of  the  upper  water  layer 
after  removing  water  molecules  within  2.8  A  of  the  indolicidin 
peptide.  The  number  of  counterions  in  the  upper  water  layer 
was  adjusted  to  maintain  charge  neutrality  in  the  presence  of 
the  positively  (+4)  charged  indolicidin. 

Calculation  of  the  Free  Energy  Profile.  For  a  detailed 
understanding  of  how  indolicidin  interacts  with  lipid  bilayers, 
we  performed  MD  simulations  with  the  peptide  z  position 
restrained  at  various  points  in  the  aqueous  phase  with  NAMD’s 
colvar  module.  We  used  a  harmonic  biasing  potential  with  a 
force  constant  of  2.0  kcal/(mol  A2)  applied  to  the  center  of 
mass  z  position  of  indolicidin.  The  z  position  bias  was 
implemented  at  every  1  A  interval  from  52  to  96  A  with  respect 
to  the  quartz  surface,  and  each  biased  simulation  lasted  25  ns. 
We  calculated  the  free  energy  or  PMF  profile  by  applying  the 
weighted  histogram  analysis  method43'44  to  the  last  15  ns  of  the 
simulation  data.  A  total  of  3.375  /is  of  the  simulation  data 
derived  from  five  different  and  independent  starting  config¬ 
urations  were  used  to  calculate  the  free  energy  profile  for  each 
bilayer  system.  However,  the  z  position  bias  for  simulations 
starting  with  the  NMR  model  structure  was  extended  up  to  125 
A  to  observe  the  structural  properties  of  the  peptide  near  the 
water/vacuum  interface. 

Lipid  Order  Parameter.  The  acyl  chain  lipid  order 
parameter  profile45  is  one  of  the  important  quantities 
describing  the  properties  of  lipid  bilayers.  To  probe  possible 
disruptions  of  the  membranes  caused  by  indolicidin,  we 
calculated  the  carbon-deuterium  lipid  order  parameter  SCD  for 


the  aliphatic  C— H  bond  vectors  of  each  carbon  atom  of  the  acyl 
chain  sn—  1  with  the  following  formula: 


SCD  ~ 


1 


(1) 


where  6  is  the  angle  between  the  C— H  vector  and  the  z  axis. 
The  order  parameter  SCD  is  zero  for  a  completely  unordered 
(isotropic)  system.  A  perfectly  ordered  acyl  chain  in  an  all-trans 
conformation  results  in  an  SCD  value  of  — O.5.45 

Electrostatic  Property  Profiles  Across  the  Interface  in 
Periodic  Boundary  Conditions.  Distribution  of  the  electric 
field  along  the  z  direction  E(z)  was  calculated  from  the  charge 
density  distribution  pq(z)  obtained  from  the  simulations  using 
the  following  relationship:28'36'46 

/-V  2Pq(z')dZ' 

e0  (2) 


E(z) 


where  e0  is  the  vacuum  permittivity  and  Lz  is  the  length  of  the 
simulation  cell  in  the  z  direction. 


■  RESULTS  AND  DISCUSSIONS 

Structural  and  Electrostatic  Properties  of  Quartz- 
Supported  Neat  Lipid  Bilayers.  To  delineate  the  underlying 
differences  between  model  eukaryotic  and  prokaryotic  cell 
membranes,  we  analyzed  structural  and  electrostatic  properties 
of  hydrated  DMPC  and  mixed  DMPC/G  bilayers  supported  on 
a  quartz  crystal  in  the  absence  of  the  added  peptide.  Figure  2 
compares  mass  density  distributions  of  lipid,  water,  and  sodium 


z(  A) 


Figure  2.  Mass  density  distributions  of  components  of  the  quartz 
crystal  supported  lipid  bilayers:  (a)  DMPC  bilayer  without  salt  ions, 
(b)  DMPC  bilayer  with  added  salt  ions,  and  (c)  DMPC/G  bilayer 
with  added  salt  ions.  Black,  red,  blue,  and  green  lines  represent  lipid, 
water,  sodium,  and  chloride,  respectively.  The  z  position  in  this  and 
subsequent  figures  is  with  respect  to  the  membrane  center. 
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and  chloride  ions  in  the  z  direction  normal  to  the  lipid  bilayer 
plane.  Density  distributions  of  lipids  and  water  (Figure  2)  were 
all  similar  regardless  of  the  compositions  of  lipids  and  presence 
or  absence  of  salt  ions.  However,  in  the  upper  bulk  water  phase 
at  z  >  0  A,  part  b  and  c  of  Figure  2  show  a  noticeable  difference 
in  sodium  ion  localization  and  penetration  into  the  model 
membranes,  with  a  deeper  sodium  ion  penetration  and 
enhanced  density  in  the  DMPC/G  lipid  bilayer  compared  to 
the  DMPC  bilayer.  Chloride  ions  exhibited  the  opposite 
behavior,  with  enhanced  density  at  the  DMPC  water  interface 
and  decreased  density  at  the  DMPC/G  water  interface.  At 
distances  farther  away  from  the  interface  (z  >  28  A),  the  nearly 
uniform  ion  distributions  in  parts  b  and  c  of  Figure  2  indicate 
the  limit  of  the  bilayer  influence  on  ion  ordering  in  the  solvent 
phase.  The  unequal  distribution  of  sodium  and  chloride  ions  at 
the  interface  resulted  in  a  slight  excess  of  a  net  system  positive 
charge  near  the  DMPC  membrane  and  a  slightly  positive 
electric  field  of  0.3  mV/A  in  the  bulk  phase  of  the  upper  water 
layer.  Part  a  of  Figure  3  shows  both  the  fluctuations  in  the  total 


z  (A) 


(b) 


z  (A) 

Figure  3.  Distributions  of  electric  field  and  water  polarization  in  the 
upper  water  layer.  Blue  lines  show  results  for  the  DMPC  bilayer 
without  salt  ions.  Black  and  red  lines  show  results  for  DMPC  and 
DMPC/G  bilayers  with  added  salt  ions,  respectively,  (a)  Electric  field 
E(z)  calculated  by  eq  2.  Dotted  lines  represent  raw  results,  which 
fluctuate  around  zero.  To  capture  the  trends  in  these  highly  fluctuating 
electric  field  values,  we  calculated  20  A  running  averages  of  E(z)  near 
the  bulklike  region  of  the  upper  water  layer  shown  as  solid  lines,  (b) 
Orientational  polarization  of  water  estimated  by  <cos0dipolez>,  the 
average  cosine  of  the  angle  between  the  water  dipole  and  the  z  axis. 


electric  field  distributions  and  the  smaller  but  consistently 
nonzero  average  electric  field  values  calculated  using  eq  2.  The 
net  positive  charge  near  the  model  eukaryotic  cell  membrane  in 
the  presence  of  salt  ions  can  play  a  role  in  the  selective 
antimicrobial  action  by  deterring  cationic  antimicrobial 
peptides  from  approaching  the  bilayer.  However,  the  magnitude 
of  the  net  effect  is  small  and  could  be  subject  to  the  specifics  of 
the  lipid  and  ions  of  the  particular  system  modeled.  Part  c  of 
Figure  2  shows  that  the  additional  sodium  counterions 
associated  with  the  anionic  lipids  are  preferentially  located 
near  the  mixed  DMPC/G  bilayer.  However,  the  electric  field 
distribution  shown  in  part  a  of  Figure  3  exhibits  a  slightly 


negative  electric  field  value  of  —0.2  mV/  A  in  the  bulk  phase  of 
the  upper  water  layer  for  the  DMPC/G  lipid  bilayer.  This 
implies  a  slightly  negative  net  effective  charge  near  the  DMPC/ 
G  lipid  bilayer  caused  by  incomplete  shielding  of  the  negative 
charges  of  anionic  lipids,  despite  the  sodium  ion's  tendency  to 
penetrate  deeper  into  the  membrane.  Importantly,  in  the 
absence  of  any  salt  ions  (part  a  of  Figure  2),  the  average  electric 
field  in  the  bulk  phase  of  water  is  nearly  zero  (part  a  of  Figure 
3).  This  confirms  that  the  presence  of  counterions  and 
physiological  salts  play  a  significant  role  in  modulating  the 
electrostatic  properties  of  the  lipid  bilayer  system. 

In  spite  of  the  fluctuations  shown  in  Figure  3,  the  electric 
fields  have  a  non-negligible  effect  in  terms  of  polarizing  the  bulk 
aqueous  solution.  To  characterize  this  effect,  we  measured 
orientational  polarization  of  water  by  calculating  the  average 
value  of  cosine  of  the  angle  between  the  water  dipole  and  the  z- 
axis,  (cos6>dipolez).  This  quantity  is  a  sensitive  probe  of  the  local 
electric  field  in  simulated  systems.28,46,47  Part  b  of  Figure  3 
shows  (cos6>dipolez)  as  a  function  of  the  z  position  for  different 
bilayer  systems.  In  the  absence  of  salt  ions,  water  was  not 
polarized,  and  (cos0dipole;Z)  was  close  to  zero  in  the  bulk  phase 
of  water  with  the  DMPC  bilayer,  consistent  with  the  net  zero 
electric  field  in  the  bulk  phase  shown  in  part  a  of  Figure  3.  On 
the  other  hand,  slightly  positive  and  negative  (cos6>dipolez)  values 
were  observed  for  the  DMPC  and  DMPC/G  bilayers, 
respectively,  consistent  with  the  electric  field  profiles  in  part  a 
of  Figure  3.  Both  the  zwitterionic  DMPC  charges  and  the 
negative  charges  introduced  by  the  anionic  DMPG  lipids  are 
strongly  shielded  by  the  counterions  and  salts  near  the  lipid 
bilayer.  However,  in  the  presence  of  salt,  the  net  small 
polarization  of  the  bulk  water  results  in  opposite  effects  in  the 
two  bilayers:  in  the  DMPC  system,  water  dipoles  are  oriented 
away  from  the  interface,  whereas,  in  the  DMPC/G  systems, 
they  are  oriented  toward  the  interface.  As  indicated  in  part  b  of 
Figure  3,  this  effect  is  small  in  magnitude  but  significant  in 
differentiating  the  electrostatic  behavior  of  the  two  systems  as  it 
relates  to  indolicidin  adsorption. 

Interaction  of  Indolicidin  with  Quartz-Supported 
Lipid  Bilayers.  To  investigate  the  interaction  between 
indolicidin  and  the  model  lipid  bilayers,  we  performed  a  series 
of  simulations  where  the  z  position  of  indolicidin  was  restrained 
at  1  A  intervals  from  the  bulk  phase  of  the  upper  water  layer  to 
the  membrane  interface.  As  described  in  the  Methods  section, 
we  used  five  different  initial  conformations  of  indolicidin  and 
two  different  compositions  of  lipid  bilayers  to  generate  results 
that  were  independent  of  starting  configuration.  Figure  4  shows 
the  resulting  PMF  profiles  combined  from  the  five  independent 
simulations.  The  placement  of  indolicidin  at  the  membrane 
surface  was  associated  with  a  minimum  PMF  of  —3.5  ±  0.2 
kcal/mol  at  z  =  25.0  A  for  the  prokaryotic  DMPC/G  bilayer, 
compared  to  —0.6  ±  0.6  kcal/mol  at  z  =  28.1  A  for  the 
eukaryotic  DMPC  bilayer.  Consistent  with  the  electric  field 
profiles  and  orientational  polarization  of  bulk  water  shown  in 
Figure  3,  indolicidin  exhibited  a  slight  attraction  to  the  mixed 
DMPC/ G  lipid  bilayer  but  was  weakly  repulsed  by  the  DMPC 
bilayer.  This  difference  in  the  PMF  profiles  of  indolicidin  with 
DMPC  and  DMPC/G  lipid  bilayers  reflects  the  selective 
association  of  indolicidin  with  the  bacterial  cell  membrane. 

Because  the  microscopically  distance-dependent  PMF  for 
indolicidin  is  not  readily  available  from  experiments,  we 
estimated  two  macroscopic  quantities  from  our  simulation 
data.  Applying  the  formalism  from  Gray  and  co-workers,48,49 
we  used  the  distance-dependent  PMF  in  Figure  4  to  estimate 
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Figure  4.  Potential  of  mean  force  (PMF)  or  free  energy  profiles  of 
indolicidin  as  a  function  of  the  distance  of  the  peptide  from  the 
membrane  center  of  DMPC  and  DMPC/G  lipid  bilayers.  PMF 
distributions  with  DMPC  and  DMPC/G  lipid  bilayers  are  shown  in 
solid  lines  with  circles  and  squares,  respectively.  Errors  reported  were 
based  on  the  95%  confidence  interval  estimated  as  1.96  times  the 
standard  error.  Standard  errors  were  estimated  from  the  analysis  of  five 
independent  PMF  distributions  obtained  from  five  sets  of  simulations 
with  different  starting  configurations  of  indolicidin  for  each  bilayer. 
The  insert  shows  the  corresponding  probability  distribution  P(z )  of 
indolicidin  on  the  two  bilayers.  The  distribution  for  the  DMPC/G 
bilayer  was  reduced  by  a  factor  of  40  for  easier  comparison  with  the 
DMPC  bilayer. 


the  free  energy  of  indolicidin  adsorption  (AG°)  to  be  —1.94  ± 
0.56  kcal/mol  for  the  model  prokaryotic  cell  membrane  and  of 
—0.05  ±  0.46  kcal/ mol  for  the  eukaryotic  membrane.  From  the 
PMF  curve  in  Figure  4,  we  further  calculated  the  maximum 
adsorption  binding  force  to  be  21  pN  at  z  =  31  A  and  41  pN  at 
z  =  30  A  for  the  DMPC  and  DMPC/ G  bilayers,  respectively. 
The  DMPC  binding  force  value  was  close  to  the  value 
determined  for  bovine  lactoferricin  adsorption  on  the  closely 
related  zwitterionic  l-palmitoyl-2-oleoyl-sn-glycero-3-phospho- 
choline  (POPC)  pure  bilayer.48 

To  characterize  the  overall  energetic  changes  as  indolicidin 
approaches  the  lipid  bilayer,  we  calculated  the  average  potential 
energy  (enthalpy)  as  a  function  of  a  center- of-mass  distance  of 
indolicidin  from  the  membrane  center.  Part  a  of  Figure  5  shows 
the  results  for  both  bilayer  systems.  The  total  potential  energy 
of  the  system  decreased,  that  is,  contributed  favorably  to  the 
free  energy  of  adsorption,  as  indolicidin  approached  both  types 
of  lipid  bilayers.  This  implies  that  the  approach  of  indolicidin 
toward  either  lipid  bilayer  was  an  enthalpically  favorable 
process.  Parts  b  and  c  of  Figure  5  show  the  further  division  of 
the  potential  energies  in  terms  of  electrostatic  and  van  der 
Waals  components,  respectively.  Counterintuitively,  the  total 
electrostatic  energies  of  the  system  did  not  depend  strongly  on 
distance  and  was  not  attractive  as  indolicidin  approached  the 
membrane;  in  fact,  the  electrostatic  component  of  the  enthalpy 
became  repulsive  at  closer  distances.  However,  the  van  der 
Waals  energies  decreased  significantly  as  indolicidin  moved 
closer  to  the  lipid  bilayers.  Part  d  of  Figure  5  shows  the 
distance-dependent  profiles  of  the  entropic  contribution  to  the 
PMF  as  estimated  by  —TAS  =  PMF  —  AD.  The  contribution  of 
the  entropic  component  (—TAS)  to  the  free  energy  of  both 
bilayer/ indolicidin  systems  was  opposite  the  contribution  of  the 
enthalpic  term.  As  indolicidin  starts  interacting  more  closely 
with  the  headgroups  of  the  lipids  and  penetrating  into  the  outer 
portion  of  the  bilayer,  there  is  an  entropic  penalty.  The  detailed 
balance  between  these  forces  determined  the  adsorption  free 
energy  and  the  locations  of  the  PMF  minima  seen  in  Figure  3. 
Although  the  statistical  uncertainty  of  the  separated  free  energy 


Figure  5.  Contribution  of  different  energy  terms  as  a  function  of 
indolicidin  z  position  from  the  membrane  center:  (a)  total  potential 
energy,  (b)  electrostatic  energy,  (c)  van  der  Waals  energy,  and  (d) 
entropic  contribution.  Solid  and  dashed  lines  represent  results  with 
DMPC  and  DMPC/ G  lipid  bilayers,  respectively.  Errors  reported  were 
based  on  95%  confidence  intervals  estimated  as  1.96  times  the 
standard  error.  Standard  errors  were  estimated  by  the  analysis  of 
energy  values  obtained  from  five  sets  of  simulations  with  different 
starting  configurations  of  indolicidin  for  each  bilayer. 


components  is  larger  than  for  the  PMF  calculation,  we 
estimated  the  AD  and  —TAS  for  both  bilayer  systems  at 
their  respective  minima  to  be  —2  kcal/mol  and  +1  kcal/mol  for 
the  DMPC  bilayer  and  —1  kcal/mol  and  —2  kcal/mol  for  the 
mixed  DMPC/G  bilayer.  Qualitatively,  the  enthalpic  changes 
are  similar,  whereas  the  entropic  terms  favor  the  DMPC/G 
system.  To  elucidate  the  nature  of  the  indolicidin  interaction 
with  the  bilayers,  we  next  characterized  the  structural  properties 
tied  to  these  energetic  changes. 

Structural  Properties  of  Indolicidin.  We  characterized 
structural  properties  of  indolicidin  as  a  function  of  its  distance 
to  the  lipid  bilayers.  Figure  6  shows  snapshots  of  indolicidin 
from  the  simulations  at  different  z  positions  with  the  mixed 
DMPC/G  lipid  bilayer.  In  part  a  of  Figure  6,  indolicidin  is 
located  in  the  middle  of  the  upper  water  layer  (z  =  48.5  A  from 
the  membrane  center).  Part  b  of  Figure  6  shows  indolicidin  at  z 
=  34.7  A  from  the  membrane  center,  where  two  adjacent 
arginine  residues  near  the  C-terminus  of  indolicidin  have  begun 
to  interact  with  the  mixed  lipid  bilayer.  Part  c  of  Figure  6  shows 
indolicidin  at  z  =  25.6  A,  where  positively  charged  lysine  and 
arginine  residues  interact  with  the  mixed  lipid  bilayer.  Part  d  of 
Figure  6  shows  indolicidin  at  z  =  19.1  A,  where  several 
hydrophobic  tryptophan  residues  are  interacting  with  the  lipid 
hydrocarbon  chains  and  several  of  the  positively  charged 
residues  are  making  contacts  with  lipid  headgroups  in  the 
interfacial  region. 
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(a)  (b)  (c)  (d) 


Figure  6.  Snapshots  of  indolicidin  at  (a)  48.5,  (b)  34.7,  (c)  25.6,  and  (d)  19.1  A  away  from  the  center  of  the  DMPC/G  lipid  bilayer.  DMPC  and 
DMPG  lipids  are  represented  by  gray  and  red  lines,  respectively.  Phosphorus  atoms  in  the  headgroups  of  DMPC  and  DMPG  lipids  are  shown  as  gray 
and  red  spheres,  respectively.  Indolicidin  is  represented  by  a  ribbon  in  cyan.  Non-hydrogen  atoms  of  positively  charged  lysine  and  arginine  residues 
and  N-terminus  are  explicitly  represented  as  balls  and  sticks.  Carbon,  nitrogen,  and  oxygen  atoms  are  shown  in  cyan,  blue,  and  red,  respectively. 
Indole  rings  of  hydrophobic  tryptophan  residues  are  represented  by  thick  gray  bonds.  Oxygen  atoms  of  water  molecules  are  shown  as  small  pink 
balls. 


Part  a  of  Figure  7  shows  the  distributions  of  backbone  rmsd 
of  indolicidin  with  respect  to  the  experimental  NMR  structure 
averaged  over  all  configurations  in  the  production  runs.  For 
both  lipid  bilayers,  the  structure  of  indolicidin  deviated 
significantly  from  the  experimental  NMR  structure.  Part  b  of 
Figure  7  shows  backbone  rmsd  of  indolicidin  as  a  function  of 
indolicidin  z  position  from  the  membrane  center.  The  average 
backbone  rmsd  values  of  indolicidin  in  bulk  water  phase  and 
near  the  bilayer  were  similar,  whereas  indolicidin  exposed  to  the 
vacuum  (z  >  60  A)  displayed  larger  average  backbone  rmsd 
values.  Part  c  of  Figure  7  shows  orientational  distributions  of 
indolicidin  measured  by  the  average  cosine  of  the  angle 
between  the  z  axis  and  the  long  axis  of  indolicidin  determined 
from  the  inertia  tensor  (<cos0peptidez>).  The  angle  0peptide>z 
ranged  from  0  to  90  degrees  with  respect  to  the  z  axis  and 
cos^peptidez  ranged  from  0  (parallel  to  the  membrane  interface) 
to  1  (perpendicular  to  the  membrane  interface).  When 
indolicidin  was  located  in  the  bulk  water  phase,  <cos0peptidez> 
was  about  0.5,  consistent  with  an  isotropic  distribution.  Near 
the  water/vacuum  interface  (60  A  <  z  <  70  A),  indolicidin 
adopted  a  perpendicular  orientation  with  respect  to  the 
interface.  More  importantly,  however,  when  indolicidin  was 
located  closer  to  the  lipid  bilayer,  <cos0peptidez>  was  slightly 
increased  (30  A  <  z  <  40  A)  but  decreased  in  region  where  the 
minima  of  the  PMF  are  located  (Figure  3).  Thus,  in  the 
energetically  most  favorable  configurations,  indolicidin  adopts  a 
roughly  parallel  orientation  with  respect  to  the  membrane 
interface.  Part  c  of  Figure  6  shows  an  example  of  such  a 


conformation.  The  structural  variations  and  orientational 
preferences  of  indolicidin  were  qualitatively  similar  for  both 
the  DMPC  and  DMPC/G  bilayer  systems. 

Khandelia  and  Kaznessis18  observed  a  persistent  cation— /r 
interaction  between  TRP11  and  ARG13  of  indolicidin  in  the 
zwitterionic  DPC  micelle.  To  determine  whether  we  observe  a 
similar  stabilizing  interaction,  we  calculated  the  distributions  of 
the  distance  between  TRP11  and  ARG13  of  indolicidin  near 
DMPC  and  DMPC/G  lipid  bilayers  as  shown  in  part  d  of 
Figure  7.  A  larger  population  of  indolicidin  with  the  TRP11- 
ARG13  distance  less  than  5  A  was  observed  near  the 
zwitterionic  DMPC  bilayer  than  the  mixed  DMPC/G  bilayer, 
which  is  consistent  with  the  simulation  results  by  Khandelia  and 
Kaznessis.18  However,  no  persistent  cation— ;r  interaction  was 
observed  for  indolicidin  near  either  lipid  bilayer  in  our 
simulations. 

Distributions  of  Positively  Charged  Residues  and 
Tryptophan  Residues.  In  Figure  8,  we  compared  distribu¬ 
tions  of  positively  charged  residues  and  hydrophobic 
tryptophan  residues  from  indolicidin  adsorbed  on  DMPC  and 
DMPC/G  lipid  bilayers.  To  select  the  indolicidin  interface- 
contact  conformations,  we  used  configurations  where  the 
center-of-mass  of  indolicidin  was  located  between  20  and  27  A 
away  from  the  membrane  center.  Part  a  of  Figure  8  shows  that 
positively  charged  residues  were  more  preferentially  distributed 
near  the  headgroups  of  the  DMPC/ G  lipid  bilayers  than  in  the 
DMPC  system.  However,  part  b  of  Figure  8  shows  that  there 
was  no  significant  difference  between  distributions  of 
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Figure  7.  Structural  properties  of  indolicidin  in  the  lipid  bilayer 
systems.  Solid  and  dashed  lines  represent  results  from  simulations  with 
DMPC  and  DMPC/G  lipid  bilayers,  respectively.  Backbone  rmsd  of 
indolicidin  was  calculated  with  respect  to  the  structure  of  indolicidin 
bound  to  a  dodecylphosphocholine  micelle  as  determined  by  NMR.* * * * * * * 8 
(a)  Distributions  of  rmsd  from  all  production  runs,  (b)  Average  rmsd 
as  a  function  of  indolicidin  z  position  from  the  membrane  center,  (c) 
Orientation  of  indolicidin  peptide  measured  by  <cos 6  peptide;z>;  the 
average  cosine  of  the  angle  between  the  peptide  long  axis  and  the 
membrane  normal,  as  a  function  of  indolicidin  z  position.  The  data  for 
z  >  60  A  corresponds  to  indolicidin  near  the  water/ vacuum  interface 
centered  around  z  =  80  A.  (d)  Distributions  of  the  distance  between 
TRP11  and  ARG13  of  indolicidin  near  DMPC  and  DMPC/G  lipid 
bilayers.  The  distance  between  TRP11  and  ARG13  (dTRP11.ARG13)  was 
defined  as  the  distance  between  the  CE2  atom  of  the  tryptophan 
indole  ring  and  the  CZ  atom  of  arginine  as  defined  in  the  PDB  file. 
The  distributions  were  calculated  with  the  conformations  of 
indolicidin  whose  z  position  is  between  20  A  and  27  A  away  from 
the  membrane  center. 


tryptophan  residues  near  the  DMPC  and  DMPC/G  lipid 
bilayers.  The  primarily  hydrophobic  interactions  stemming 

from  tryptophan/lipid  contacts  were  roughly  similar  in  nature, 

commensurate  with  the  similar  enthalpic  profiles  shown  in 

Figure  5  for  both  types  of  bilayers.  Thus,  the  distinction 
between  the  two  bilayers  in  terms  of  interactions  lies  primarily 

in  how  the  positively  charged  residues  are  distributed  in  the 
interface  region  of  the  bilayer  and  how  the  system  of  solvent 

and  ion  molecules  accommodate  these  configurations.  Given 
that  the  enthalpic  component  derived  from  electrostatic 
interactions  between  the  two  systems  is  not  that  different  due 
to  charge  conservation  and  counterion  shielding,  the 

conformations  underlying  both  distributions  in  part  a  of  Figure 

8  are  enthalpically  equivalent.  The  difference  in  density  of 
positively  charged  residues  located  in  the  bilayer  headgroup 
region  around  z  =  20  A  in  part  a  of  Figure  8  instead  points  to  a 
mechanism  whereby  positively  charged  ions  normally  asso¬ 
ciated  with  the  bilayer  are  replaced  by  positive  charges  from 
indolicidin.  This  type  of  an  effect  would  essentially  be 
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Figure  8.  Distributions  of  positively  charged  residues  and  hydrophobic 
tryptophan  residues.  Distributions  were  calculated  from  configurations 
where  indolicidin  was  located  between  20  and  27  A  away  from  the 
membrane  center,  (a)  Distributions  of  positively  charged  residues  were 
estimated  as  the  sum  of  density  distributions  of  the  N  atom  at  the  N- 
terminus,  the  NZ  atom  of  lysine,  and  the  CZ  atom  of  arginine  as 
defined  in  the  PDB  file,  (b)  Distributions  of  tryptophan  residues 
calculated  by  the  density  distribution  of  the  CE2  atom  of  the  indole 
ring  in  the  PDB  file. 


enthalpically  neutral  but  entropically  favorable.  This  interpre¬ 
tation  is  consistent  with  the  observed  difference  in  entropy 
changes  seen  for  the  two  different  bilayers. 

These  results  indicate  that  the  electrostatic  interactions 
between  positively  charged  residues  of  indolicidin  and  the 
headgroups  of  anionic  lipids  in  the  mixed  lipid  bilayer  underlie 
the  specificity,  but  not  the  strength  of  attraction,  of  indolicidin 
toward  the  bacterial  cell  membrane.  Instead,  the  increased 
strength  of  adsorption  free  energy  in  the  DMPC/G  model 
prokaryotic  membrane  versus  the  DMPC  eukaryotic  membrane 
is  postulated  to  be  derived  from  a  gain  in  entropy  due  to 
counterion  release.49,50  Figure  9  shows  the  positive  charge 
distribution  in  both  bilayer  systems,  both  in  the  presence  and 
absence  of  indolicidin  at  the  bilayer  surface.  Part  a  of  Figure  9 
shows  that  the  total  charge  distribution  seen  in  the  neat  system 
has  been  replaced  by  a  combination  of  indolicidin  charges  and 
Na+  atoms  for  the  DMPC/G  bilayer.  Consequently,  sodium 
atoms  strongly  localized  to  the  surface  in  the  absence  of 
indolicidin  were  released  into  the  bulk  solution  upon  surface 
adsorption  of  indolicidin,  providing  a  relative  gain  in  entropy. 
For  the  DMPC  bilayer  however,  parts  b  and  c  of  Figure  2  show 
that  the  sodium  atoms  in  the  neat  solution  are  not  strongly 
localized  to  the  bilayer  surface  compared  to  the  bulk  solution. 
Consequently,  although  adsorption  of  indolicidin  at  the  surface 
also  released  sodium  atoms,  it  did  not  provide  the  same  relative 
gain  in  entropy  as  for  the  DMPC/G  bilayer.  The  observed 
entropy  change  is  an  intrinsic  property  of  the  entire  system  and 
cannot  strictly  be  decomposed  into  atomic  components. 
However,  in  support  of  our  interpretation,  Vivcharuk  and 
Kaznessis50  found  the  peptide— membrane  attraction  to  be 
dominated  by  the  entropy  increase  due  to  the  release  of 
counterions  in  a  POPG/POPE  lipid  bilayer  with  MD 
simulations  where  orientations  and  conformations  of  the 
peptide  were  constrained. 
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Figure  9.  Distributions  of  Na+  and  positively  charged  residues  in  the 
presence  and  absence  of  indolicidin.  Distributions  for  the  systems  with 
indolicidin  adsorbed  onto  the  bilayer  were  calculated  from 
configurations  where  the  peptide  was  located  between  20  and  27  A 
away  from  the  membrane  center.  The  contribution  from  the  positive 
charged  residues  of  indolicidin  was  calculated  as  the  sum  of  density 
distributions  of  the  N  atom  at  the  N-terminus,  the  NZ  atom  of  lysine, 
and  the  CZ  atom  of  arginine  as  defined  in  the  PDB  file,  (a)  The  total 
positive  charge  distribution  for  the  DMPC/G  bilayer  including  Na+ 
and  indolicidin  charges  (He,  Lys,  Arg,  Na+  DMPC/G/IL),  the  charge 
distribution  due  to  Na+  only  in  the  presence  of  indolicidin  (Na+ 
DMPC/G/IL),  and  the  charge  distribution  due  to  Na+  in  the  neat 
systems  taken  from  Figure  2  (Na+  neat  DMPC/G).  (b)  The  same 
distributions  for  the  DMPC  bilayer. 

Bilayer  Perturbations  and  Lipid  Order  Parameters. 

Lastly,  we  examined  the  lipid  bilayer  itself  for  structural  features 
that  could  differentiate  the  response  of  the  bilayers  to 
indolicidin  adsorption.  We  created  profiles  of  the  lipid  order 
parameter  —  SCD  along  the  lipid  acyl  chain  for  DMPC  and 
DMPC/G  bilayers  using  eq  1.  We  constructed  —  SCD  profiles 
for  phospholipids  that  were  in  contact  with  indolicidin  and 
compared  them  with  —  SCD  profiles  obtained  when  the  peptide 
was  sufficiently  distant  (z  >  50  A  from  the  membrane  center) 
such  that  no  peptide-phospholipid  contact  existed.  A  lipid 
molecule  in  the  upper  leaflet  of  the  bilayer  was  considered  in 
contact  with  the  peptide  if  the  center-of-mass  z  position  of 
indolicidin  was  between  20  and  27  A  away  from  the  membrane 
center  and  the  radial  distance  of  the  phosphorus  atom  of  the 
lipid  from  the  cylindrical  axis  along  the  z  direction  that  passes 
through  the  center-of-mass  position  of  peptide  was  within  8  A. 
Parts  a  and  b  of  Figure  10  show  the  calculated  —  SCD  profiles  for 
the  DMPC  and  DMPC/G  bilayers,  respectively.  For  — 1 Scd 
profiles  without  the  peptide— lipid  contact,  the  profiles 
combined  from  both  upper  and  lower  leaflets  of  lipid  bilayers 
are  shown.  We  did  not  observe  any  significant  difference 
between  —  SCD  profiles  calculated  separately  for  upper  and 
lower  leaflets  of  bilayers  (data  not  shown).  This  indicates  that 
the  effect  of  the  quartz  solid  support  on  the  structural 
properties  of  the  lipid  bilayer  was  negligible  when  there  was  a 
thick  water  layer  between  the  quartz  surface  and  the  lipid 
bilayer  as  we  observed  in  our  previous  study.28  As  expected, 
— SCD  values  from  lipids  in  contact  with  the  peptide  were 
smaller  than  lipids  not  in  contact,  indicating  a  more  disordered 
acyl  chain  for  lipids  in  contact  with  indolicidin.  The  average  z 
displacements  of  the  phosphorus  atoms  of  lipids  in  contact  with 
the  peptide  from  those  of  upper-leaflet  lipids  not  in  contact 


Figure  10.  Order  parameter  —  SCD  along  the  lipid  acyl  chain  of  (a) 
DMPC  and  (b)  DMPC/G  bilayers  calculated  with  eq  1.  The  carbon 
atoms  are  numbered  sequentially  from  head  to  tail  regions  of  a  lipid. 
Order  parameters  for  lipids  with  and  without  the  peptide— lipid 
contact  are  represented  by  lines  with  and  without  circles,  respectively. 

with  the  peptide  were  —1.69  and  —2.0  A  for  DMPC  and 
DMPC/G  bilayers,  respectively.  Even  though  —  SCD  values  for 
carbon  numbers  3  through  9  near  the  lipid  headgroup  in  the 
DMPC/G  bilayer  decreased  slightly  more  than  those  for  the 
DMPC  bilayer  in  contact  with  the  peptide,  we  found  no 
significant  difference  in  changes  in  —  SCD  profiles  between 
DMPC  and  DMPC/G  lipid  bilayers. 

■  SUMMARY  AND  CONCLUSIONS 

In  this  study,  we  characterized  the  interaction  of  the 
antimicrobial  peptide  indolicidin  with  solid- supported  model 
bacterial  and  eukaryotic  cell  membranes.  We  based  our  study 
on  extensive  all-atom  molecular  dynamics  simulations  com¬ 
bined  with  rigorous  treatment  of  the  electrostatic  boundary 
condition  to  derive  structural  and  energetic  insights.  The 
calculated  free  energy  profiles  of  indolicidin  in  proximity  to 
model  bacterial  and  eukaryotic  cell  membranes  revealed  that 
there  was  a  preferential  attraction  to  the  bacterial  cell 
membrane.  Placement  of  indolicidin  at  the  membrane  surface 
was  associated  with  a  minimum  PMF  of  —3.5  ±  0.2  kcal/mol  at 
the  prokaryotic  membrane,  compared  to  —0.6  ±  0.6  kcal/mol 
at  the  eukaryotic  membrane.  We  calculated  the  free  energy  of 
adsorption  (AG°)  from  the  PMF  distance  profile  to  be  —1.94 
±  0.56  kcal/ mol  for  the  model  prokaryotic  cell  membrane  and 
—0.05  ±  0.46  kcal/mol  for  the  model  eukaryotic  membrane. 
We  also  observed  that  indolicidin  was  preferentially  located 
near  the  membrane  interface  lying  parallel  to  the  interface  in 
agreement  with  previous  experimental  and  computational 
studies. 

We  analyzed  the  energetic  components  and  the  structural 
changes  for  both  types  of  bilayers  to  gain  insights  into  the 
nature  of  indolicidin  adsorption.  Although  we  cannot  fully 
address  all  possible  lytic  mechanisms  indolicidin  is  engaged  in, 
the  initial  association  of  the  peptide  with  prokaryotic 
membranes  is  a  required  step  for  biological  activity.  We 
found  that  the  effect  of  counterions  and  salts  at  physiological 
concentrations  minimized  the  direct  electrostatic  component  of 
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the  enthalpic  contributions  to  the  PMF.  Instead,  indolicidin 
derived  the  bulk  attractive  force  from  van  der  Waals 
interactions.  This  attractive  force  was  counteracted  by  entropic 
penalties  as  the  indolicidin  molecule  approached  the  bilayer 
interface  from  the  solution.  These  effects  were  roughly  similar 
for  both  model  membranes  systems,  but  with  one  important 
difference.  The  preferential  placement  of  the  positively  charged 
residues  close  to  lipid  headgroups  of  the  DMPC/G  bilayer  as 
opposed  to  the  DMPC  bilayer  allowed  more  counterions  to  be 
released  into  the  bulk  aqueous  solution,  resulting  in  a  reduced 
entropic  penalty. 

In  the  current  study,  we  investigated  adsorption  of  a  single 
indolicidin  molecule  on  to  the  bilayer.  We  fixed  the  simulation 
box  sizes  in  directions  parallel  to  the  membrane/water  interface 
to  maintain  the  geometry  of  the  quartz  (Oil)  surface.  For  a 
broader  understanding  of  membrane  perturbations  at  deeper 
penetrations  and  higher  concentrations51  of  indolicidin,  as  well 
as  other  antimicrobial  peptide  classes,  we  would  require  a 
flexible  simulation  cell  that  can  accommodate  larger  rearrange¬ 
ment  of  the  bilayer  than  seen  here.  We  are  currently  developing 
a  constant  pressure  simulation  methodology  to  accommodate 
larger  lipid  bilayer  perturbations  due  to  antimicrobial  peptides. 

In  summary,  the  nature  of  the  indolicidin  surface  adsorption 
depended  on  an  electrostatic  guiding  component,  an  attractive 
enthalpic  component  derived  from  van  der  Waals  interactions, 
and  an  entropic  factor  dependent  on  the  amount  of  counterion 
release  from  the  bilayer  surface.  Commensurate  with  this  view, 
introducing  additional  charged  residues  into  indolicidin  will  not 
greatly  alter  the  strength  of  adsorption  but  could  alter  the  type 
of  bilayer  the  peptide  would  adsorb  on.17  It  also  follows  that 
modulating  the  overall  van  der  Waals  interactions  by 
substituting  the  bulkier  tryptophan  residues  with  smaller 
phenylalanine  residues  will,  however,  noticeably  decrease  the 
adsorption  free  energy.17  Engineering  indolicidin-based  or 
other  cationic  peptides  based  on  maintaining  or  improving 
the  differential  surface  adsorption  profile  between  eukaryotic 
and  prokaryotic  membranes  should  utilize  these  principles, 
whereas  the  specifics  of  improving  bacterial  lysis  and  reducing 
hemolytic  activity  could  be  guided  by  different  principles. 
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